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Abstract. This contribution presents a derivation of the steady-state distribution of velocities and dis- 
tances of driven particles on a onedimensional periodic ring. We will compare two different situations: (i) 
symmetrical interaction forces fulfilling Newton's law of "actio = reactio" and (ii) asymmetric, forwardly 
directed interactions as, for example in vehicular traffic. Surprisingly, the steady-state velocity and dis- 
tance distributions for asymmetric interactions and driving terms agree with the equilibrium distributions 
of classical many-particle systems with symmetrical interactions, if the system is large enough. This an- 
alytical result is confirmed by computer simulations and establishes the possibility of approximating the 
steady state statistics in driven many-particle systems by Hamiltonian systems. Our finding is also useful to 
understand the various departure time distributions of queueing systems as a possible effect of interactions 
among the elements in the respective queue [D. Helbing et al., Physica A 363, 62 (2006)]. 
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1 Introduction 

Classical many-particle systems such as ideal gases are 
characterized by the applicability of Newton's laws of me- 
chanics, which particularly includes the law of "actio = 
reactio" . With these basic laws, many fundamental prop- 
erties can be derived, like the conservation of momen- 
tum and energy. Such many-particle systems are known 
as Hamiltonian systems. Even statistical physics and ther- 
modynamics are based on these relationships. 

But what would happen if the law of "actio = reac- 
■ tio" would not hold and the particle interactions would 
not fulfil momentum and energy conservation? One ex- 
ample of such a system are driven Brownian particles 
where many results for stationary distributions are avail- 
able. However, in this system, interactions are taken into 
account only implicitely by a nonlinear friction function 
[1], or by a symmetric interaction potential [2|. In other 
systems such as vehicular traffic one would like to model 
the non-symmetric interactions between the particles (ve- 
hicles) explicitely. Would it still be possible to find (ana- 
lytical) formulas for the stationary distributions of veloc- 
ities and distances? In fact, although a statistical physics 
formalism for driven systems is needed, there are still not 
many results available. The existing results mainly con- 
cern the study of traffic-jam related condensation phe- 
nomena by means of the master equation [3] or the Fokker- 
Planck equation 0]. 



These considerations have assumed certain arrival and 
departure rates to or from any forming vehicle clusters, 
but they have not explicitly represented the acceleration 
or deceleration dynamics of interacting vehicles in space 
and time, which may lead to dynamic self-organization 
phenomena such as emergent stop-and-go waves [5]. In 
the following, we will take this dynamics of interacting 
particles into account. 

In fact, we are seeking for a method to treat dissipative 
driven many-particle systems in a similar way as Hamilto- 
nian systems. The idea is that the dissipation in the sys- 
tem would be balanced by the effect of the driving force, 
at least in a closed (circular) system in the limit of large 
particle numbers. This idea has been used to evaluate the 
vehicle interaction potential [61718] . but questioned to be 
applicable to systems with asymmetric, forwardly directed 
interactions. Moreover, the method has been restricted to 
a very limited number of potentials U(s), as the normal- 
ization factor of the distance distribution could not be 
analytically determined in general. 

For onedimensional classical Hamiltonian gases, many 
results have been previously derived in the framework of 
Random Matrix Theory [9110] . According to this, if cou- 
pled to a thermal bath, the velocity distribution of gas 
particles is Gaussian and the distance distribution g(s) 
can be written as 
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where U(s) is the interaction potential, is the velocity 
variance (i.e. proportional to the temperature), and B de- 
pends on the particle density. It would be very desireable 
to have a similar result for driven many-particle systems, 
as this would allow one to determine the interaction po- 
tential and interaction force among driven particles in the 
presence of fluctuations. Our hope is that, in the station- 
ary state, the dissipative interactions and the driving term 
would somehow cancel out on average, so that the be- 
havior would approximately correspond to a Hamiltonian 
system, as assumed in the Refs. |7I8I10I11I12) . 

In fact, in this paper we will show that this idea is cor- 
rect for onedimensional systems in the limit of large par- 
ticle numbers. Even forwardly directed dissipative driven 
many-particle systems behave approximately like Hamil- 
tonian systems if they are far away from a dynamic insta- 
bility. 

Our paper is structured as follows: In the next section, 
we lay out the theoretical basis and derive the form of 
the onedimensional Hamiltonian as well as the conditions, 
under which it provides a correct description. In Sec. [3] we 
formulate the predictions in a form that can be tested by 
simulating representative many-particle systems. The ac- 
tual simulations and their results are presented in Section 
IH after which we conclude with a discussion. 



2.1 Fokker-Planck equation for velocities and distances 

In order to determine the statistical distributions of the 
velocities and distances of n particles i, it is helpful to 
rewrite the above stochastic differential equation (Langevin 
equation) in terms of an equivalent Fokker-Planck equa- 
tion. With the definitions 

s = (si, . . . ,s n ) , v = (V!, ...,v n ), (5) 
W(si, s<_i) = v + r[/(sj) - jf(st-i)} , 

and P = P(si, ...,s n ,vi,...,v n ,t) = P(s,v,t) (6) 

this Fokker-Planck equation reads [13] 

—dsi I dt 

V V ' 

=dVi/dt—£i 

where we assume periodic boundary conditions Vk+ n (t) = 
Vk(t) and Sk+n(t) — Sfe(t) for a onedimensional ring of 
length L with n particles on it. In the following, we will 
show that the ansatz 



2 Driven many-particle model with dissipative 
interactions 

In the onedimensional driven-many particle system we dis- 
cuss, point-like particles i change their location Xi(t) in 
time t according to the equation of motion 



dxi 
~~dt 



Vi(t) , 



and their temporal velocity change dvi/dt is assumed to 
be given by the following stochastic acceleration equation: 



dvj _ vp - Vj 
dt ~ t 



+ /(*) -7/(*-i) +&(*)• (3) 



Here, v denotes the "free" or "desired" velocity and £j(t) 
represents a white noise fluctuation term satisfying 



&(*)) = o 



D5 i3 S{t 



•0. 



(4) 



where D is a velocity-diffusion constant. The particle mass 
m-i has been set to 1, and / (sj) < describes a repulsive 
interaction force, which depends on the particle distance 
s,(i) = Xi+\{t) — Xi{t). The term 7/(sj_i) with < 7 < 1 
allows to study different cases: 7 = 1 corresponds to the 
classical case of symmetrical interactions in forward and 
backward direction, fulfilling the physical law of "actio = 
reactio" .7 = corresponds to the case of forwardly di- 
rected interactions only, which is, for example, applicable 
to vehicles. 



P(s,v) = P(s u . 



,v n ) 



26 



(8) 



is a stationary solution of the above Fokker-Planck equa- 
tion, if the parameters V and 9 are properly chosen, and 
if the so-called interaction potential U is defined by 



U(Si) = ^r 1 I ds f(s) 



In Eq. 



J d n s J d n v P(s,v) 



(9) 



(10) 



is the normalization constant, and the Lagrange parame- 
ter B is required to meed the constraint J7 Sj = L deter- 
mining the actual particle density. Moreover, 

V(t) - { Vi ) - Jd n S J <rvViP(B,V,t) (11) 

is the average velocity and 
6(t) = ((vi-V) 2 ) = [ d n s [d n v ( Vl -V) 2 P(s,v,t) (12) 



the velocity variance. In the following, we will restrict our 
investigation to the stationary case with dV/dt — and 
d6/dt = 0, which presupposes that the deterministic part 
of Eq. fulfils the linear stability condition 



i)2 df{L/n) ^1 + 7 



ds 



2r 2 



(13) 
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(see Ref. [14] for the method to determine this formula). 
Otherwise, dynamic patterns such as stop-and-go waves 
may emerge from the dissipative interactions of driven 
particles !5 a . Notice that in the Hamiltonian case, 7 = 1, 
the stability condition is always satisfied. Furthermore, the 
factorization assumption (|SJ) requires that all variables Si 
and Vi are statistically independent from each other. Ac- 
cording to numerical simulations, this is only the case if 
(1 — "i) 2 df /ds is much smaller than the right-hand side of 
(I13p . i.e., the system is far from the instability point. 

With the ansatz ([5]), the three terms of the Fokker- 
Planck equation J7]) can be written as 



i 



B 



1 dUjsj) 
9 dsi 

(l + 7)/(*) 
29 



P 



B 



P. 



(14) 



Y — 



d (W{s i ,s i - 1 )-Vi 



P 



(Vi - V) 



P, (15) 



and 



+ 



Vi-V 



2 dm 2 2 
We will now use the fact that 



p. 



(16) 



(17) 



for any i-dependent variable gi, i.e. indices can be shifted 
because of the assumed periodic boundary conditions. In 
this way we find 



dP _ 1 

~dt ~ 9 



E 



—^—{v i+ x - v l )f{s i )P + n \—-^Q 



-y 

9 



vq - Vi 



+ f(si) - 7/(«i-i) 



(vi - V)P 



29 



(18) 



Remarkably, this equation does not depend on the La- 
grange parameter B anymore, which is needed to adjust 
to the particle density. 

Note that ansatz J5J) can only be a stationary solution 
with dP/dt = 0, if 

This relationship corresponds to the fluctuation-dissipation 
theorem. Applying it also to the last term of Eq. lfT8|) and 



using the decompositions (i>i+i 
and (v -Vi) = (v -V)- («j 



■Vi) = (v i+1 -V)-(vi-V) 
- V), we find 



8P 

~dt ~ ^ 



ft ^ 



1-7 
29 



(v i+1 +Vi-2V)f( Si )P 
- V)(vi - V) 



-P, 



and, with the factorization assumption ((? 
indices again according to (fl"7]) , we obtain 



dP 
~dt 



1 



-E 



(20) 
| and shifting 

V) P. (21) 



We will distinguish the following cases: 

1. In the case of a classical many-particle system with 
momentum conservation (7 = 1) and energy conserva- 
tion, i.e. no dissipation (r — > 00), we find dP/dt = 0, 
i.e. ansatz (jHJ) is an exact stationary solution of the 
Fokker-Planck equation ([7]). 

2. In the case of forwardly directed interactions as in ve- 
hicle traffic (7 = 0) and vanishing correlations, we have 

lim — ^^(vj — V) 



lim -VVw* 

n— >oo 77, ' * 



lim — 



f(Si) + 



vo-V 



(22) 



The first factor vanishes because of V — lirn„_ (X > i ^ i^, 
but the second factor disappears as well: Dividing Eq. ^ 
by n and summing up over i gives 

n ^ n ^— ^ t 



71 ' 



(23) 



In the limit n — > 00 of large enough particle num- 
bers n, the left-hand side converges to dV/dt, while 
the last term on the right-hand side converges to 0. In 
the assumed stationary case with dV/ dt = and using 
vq — Vi = (vq — V) — (vi — V), this implies 







lim 

n — >oo 77, 



= lim 

n— *oo n 



n ^ 



Vq - Vi 



v -V 



•/(*) 



(24) 



because of — V . v.- = V = — Y\- V. 

We conclude that the factorisation ansatz ([8]) satisfies the 
Fokker-Planck equation J7]) if either the momentum is con- 
served (7 = 1), or if the single-particle gaps and velocities 
are independent from each other. Note that, in order to 
arrive at this conclusion, the special form §§§ for the in- 
teraction potential, particularly the prefactor (1 +7)/% is 
required. 
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2.2 Hamiltonian description 

An alternative approach is the Hamiltonian description. 
For this purpose, let us investigate the Hamiltonian 



H = T + V = J2 



(Vi - Vf 



(25) 



If dV/dt — 0, we can derive the following relations: 



dH _ dT dV 
dt dt dt 



dU( Si ) 
dsi 

dsi dx i+ i 



dxi dt dxi+i dt 
dvi x-^ 1+7 



E(^ ^) 



dt ^ 2 

i 

v - Vi 



f(si){v i+1 - Vi) 



+ f(si)-7f(si-i)+Zi(t) 



E 



1+7 



E 



1-7 



(»< + Uj+i - 2V)f{ Si ) 



E 
E 



Oo - vo(«i - v) 



(Vi - V) 



•E(«i-^ 4 (t). 



Comparing this with (|2"D|) shows that 



<9P P 1 



9* 



dt 6 



E 



(«< - vf 



(vi-VMt) 



(26) 



P. 



(27) 

Correspondingly, in the stationary state dP/dt = we 
have 



dH 
~dt 



E 



(Vi - Vf 



E^-^l^*) 



T 



(28) 



T 



F) 2 



-<[/(*i)-7/(*i-i)](«i-n) 
1 d0 vq-V 

2~db ~ 



<(vi)-V) + - 



T T 

-</(ai)-7/(«i-i))((«i)-V). (29) 

This can be found by multiplication of Eq. ([3]) with 
(«i — V) and calculation of the ensemble average, us- 
ing the factorization ansatz (JSj> . The first term on the 
right-hand side vanishes under the assumption of a sta- 
tionary state. The second and the fourth term vanish 
because of (vi) — V. Therefore, 



(ti{vi-V)) = ~, 

T 

and, together with Eq. ([28]) . we arrive at 



(fK("*-? ( 



Vi - Vf 



(30) 



0. (31) 



That is, in the statistical average we have dTL/dt = 0. 
The same is expected for the average Hamiltonian per 
particle, Pi = TL/n of systems with many particles. 
In fact, simulations show that Pi fluctuates with am- 
plitudes oc l/\/n, while the Hamiltonian Ti itself fluc- 
tates with amplitudes oc y/n, which is consistent with 
equilibrium hydrodynamic systems. As a consequence, 
stationary driven dissipative systems behave approx- 
imately Hamiltonian, even if the interactions are for- 
wardly directed and Newton's law "actio = reactio" is 
violated. This is, why the Hamiltonian statistics 



P(s,v) =JVc- H ^ e 



(32) 



In order to investigate under which conditions the Hamil- 
tonian is conserved in the statistical average, we will cal- 
culate (tjt) for two different cases: 

1. In a conservative system with no fluctuations (£j(t) = 
= D) and no dissipation (r — + oo), we have dTi/dt = 
0, independently of whether the interactions are sym- 
metric or forwardly directed. 

2. For many-particle systems with fluctuation terms and/or 
dissipation, one can show 



v -V 



(vi - V) 



(the canonical distribution) is an approximate station- 
ary solution of our driven dissipative many-particle 
system. Note that the contribution J7 Bsi = BL in 
Eq. ijHJ gives just a constant prefactor and can be ab- 
sorbed into the normalization factor. 

In conclusion, the equilibrium solution 1(5)) of conservative 
many-particle systems is also a good approximation for 
the steady-state solutions (dP/dt = 0) of driven many- 
particle systems of kind ([3]) with asymmetrical interac- 
tions, driving and dissipation effects, if the system is large 
enough, i.e. n>l (for small systems, we expect that fluc- 
tuations become essential), and if the correlations between 
the gaps and velocities are insignificant. 



3 Application to stochastic traffic models 

In this section, we will formulate the results of the previous 
section in a way that can be tested by means of simulating 
specific models. 
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3.1 Single-particle distributions 

The factorisation ((HJ can be written in the form 

n n 

P(si, ...,s n ,vi,...,v n ) = IJs(si) Y[ K v i) > ( 33 ) 

»=i j=i 

i.e., the statistics of the particles can be described by the 
single-particle gap distribution function 



g(s) = Aa-M'V e+B '\ 
and the single-particle velocity distribution 



h(v) = 



-(v-V) 2 /(28) 



(34) 



(35) 



Here, A is a normalization constant, B a Lagrangian pa- 
rameter ensuring the density constraint, V the average 
velocity, and 9 the velocity variance. With the exception 
of all quantities are dependent on the particle density 
P- 



3.2 Gap distribution 

The two constants A and B of the gap distribution ((34 
are determined using the normalisation condition 



dsg(s) = 1 



(36) 



and the constraint that the average gap is equal to the 
inverse of the global density, 



oo 

J ds sg{s) = - 



(37) 



Defining the integrals 



oo 

Io(B) = J ds exp 



(U{ t 



U{s) 



I\ (B) — J ds s exp 
o 

oo 

A,;//, j ds s 2 exp f { ( ! " i 



Bs 



Bs 



Bs 



we get 



.4 



Io(B) ' 

For B, we find the transcendental equation 

h(B) _ 1 



Ah{B) 



Io(B) P 



(38) 



(39) 



Using Newton's method with the initial guess Bq = p + 
1 ja s with a s defined in Eq. (|4"3")l , one obtains for the fc-th 
iteration 

B k+1 - B k + —— -j-, (4U) 

p{hh - It) 

where the integrals on the right-hand side are evaluated at 
B = Bk- It turns out that this method converges within 
very few iterations, unless U(s e ) = 17(1/ p) 3> 9. In this 
case, however, the second derivative U"(s c ) of the effective 
potential (J9j) typically satisfies the condition 



\U"(s c )\ (l + 7 )|/'( Se )| >>p2 ^ 



9 



Dt 



(41) 



allowing an asympotic expansion of J34[) that eventually 
leads to a Gaussian gap distribution, 



i 



C/"( Sc )|»p 2 



2~KO s 



with 



Dt 



(l+ 7 )/'(s e )- 



(42) 



(43) 



Remarkably, the ranges of applicability of (|40|) and ((42 
generally overlap, allowing a fast and robust solution. 



3.3 Velocity distribution and kinetic energy 

Equation (j3"5() states that, regardless of the density, of the 
potential, and of the directions of the interactions, the 
single-particle velocity distribution is Gaussian. The ex- 
pectation value is equal to that of the stationary velocity 
without fluctating terms. Furthermore, the velocity vari- 
ance satisfies the fluctuation-dissipation theorem (flTI)) . i.e., 
the average kinetic energy per particle should be given by 



(44) 



4 Results 

In this section, we will show by means of computer simu- 
lations, that the main predictions (|3"4"1) . ([33)l . and (|4"4"1) are 
valid if the system is in a regime that is far from any col- 
lective instability. In order to quantify this condition, we 
observe from (| 13[) that the system becomes linearly unsta- 
ble if the relaxation time r exceeds some critical value r c . 
Thus, t controls the stability properties which allows us 
to define the dimensionless reduced control parameter 



r 



(45) 



Note that r = denotes equilibrium, while the linear 
threshold is characterized by r = 1. In particular, in the 
momentum-conserving case 7 = 1, we have always r = 0. 
The condition "far away from the instability point" can 
be formulated by the condition r<l. 
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Parameter 


Value 


Desired velocity vo 


30 m/s 


Velocity relaxation time r 


0.2 s 


Interaction length k n t 


20 m 


Shape parameter j3 


0.5 


Symmetry parameter 7 


and 1 


Fluctuating force D 


20 m 2 /s 3 


Density p 


12 /km and 30 /km 



Table 1. Parameter values used in the simulations for the 
stochastic OVM © with (gB]) and (Ti7|) . 



Parameter 


Value 


Desired velocity vo 


30 m/s 


Velocity relaxation time r 


2s 


Interaction distance ii n t 


20 m 


Acceleration ao 


2 m/s 2 


Interaction exponent S 


2 


Symmetry parameter 7 


and 1 


Fluctuating force D 


0.2m 2 /s 3 


Density p 


10 /km 



Table 2. Parameter values used in the simulations for the 
stochastic power- law model <(3j with (|52[) . 



4.1 Selected models 

In order to obtain a specific model, the interaction force 
of the stochastic differential equation ([3]) has to be speci- 
fied. We will simulate two types of interaction forces that 
are based on (i) the optimal-velocity model, and (ii) on a 
power law. 



4.1.1 Stochastic optimal-velocity model (sOVM) 

In the stochastic optimal-velocity model (sOVM), the in- 
teraction force / is given by [5] 



sOVM 



Vovm(s) - Vq 



(46) 



where we assume the optimal- velocity function 

(it-P) +*anh(/3) 



Vovm(s) 



tanh i 



1 + tanh(/3) 



(47) 



Table [T] summarizes the meaning of the model parameters 
and the values used in the simulations. Notice that the 
conventional optimal- velocity model of Bando et. al. [T5] 
is obtained for the special case 7 = and D = in the 
Eqs. ((3]) and (g]), respectively. 



The sOVM has the following properties: The expecta- 
tion value of the velocity for stationary conditions is given 

by 

Kovm(s) = w + (1 - i)V vm{s). (48) 
Furthermore, the effective potential © can be calculated 
analytically, resulting in 



(49) 



U s ovm(s) = t/ m <U 4- exp 



-2(7^-/3 

tint 



where the prefactor is given by 

(1 + 7)w /int 



(50) 



2t(1 + tanh/3)' 

The dynamics of this model becomes linearly unstable if 
the relaxation time exceeds the critical value r c given by 

1 + 7 



T-c(p) 



2(1 



(51) 



7)% V m(1/p)' 
see Eq. (fT3"|). For the parameters specified in Table [TJ 7 = 
0, and p = 30/km, the critical value is given by r c (-f — 
0) = \/{2Vq YM {\/ p)) = 1.51s. This corresponds to r = 
0.132 when the parameters of Table[T]are assumed. 



4.1.2 Stochastic power law model (sPLM) 

An alternative, more physics-oriented model assumes that 
the interaction forces obey a power law [6|7j : 

/ s plm(s) = -ao ( — J , (52) 



which, together with (J2]) and ([3]), results in the stochastic 
power-law model (sPLM) . The associated effective poten- 
tial §§§ is given by 



sPLM 



(s) 



(1 + 7 W? B 



(53) 



(s - ly- 1 ' 

and the expectation value for the velocity is equal to 

KsPLm(s) = v - (1 - 7)ra (-^J ■ ( 54 ) 

This model differs qualitatively from the stochastic OVM 
in the following aspects: 

— If 7 < 1, the stationary velocity becomes zero for a 
finite average gap, s c (V — 0) = [(1 — ^rao/vo] 1 ^ ■ In 
contrast, the stationary velocity of the sOVM for any 
7 > is nonzero, even at maximum density, s = 0. 

— The potential (|53|) of the stochastic power-law model 
diverges for s — > 0, while the potential (|4"9]) of the 
sOVM remains finite. For the chosen parameters and 
7 = 0, we have IWm(0) = U \n(2) = 1347 m 2 /s 2 . 
Consequently, any car approaching a standing vehi- 
cle with a velocity exceeding ^/2[/ s ovm(0) = 52 m/s 
will lead to a rear-end collision. This velocity decreases 
with increasing values of r, reaching 18.9 m/s at the 
limit of linear stability. In contrast, no such collisions 
are possible in the stochastic power-law model. Nev- 
ertheless, this model can become linearly unstable as 
well. 
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4.2 Simulations 

We have simulated a closed ring road of length L = 9 km 
for the sOVM, and L = 40 km for the sPLM. We also have 
simulated larger systems resulting in no significant differ- 
ences. We have started the simulations with deterministic 
initial conditions Si — 1/p, and Vi = V(si), corresponding 
to a single-particle distribution function 



Pi(s,v,0) = 6 (v-v(- 



d | a- - 

P 



(55) 



where S(.) represents Dirac's delta function. Since this ini- 
tial condition does not correspond to a stationary solution, 
we have run the simulations for a transient time of 72 000 s, 
before recording the results for further 36 000 s. 

For the numerical update, we have applied the explicit 
scheme 



Vi(t + At) = Vi{t) + cii(t)At + z t vDAt, (56) 

~vi(t) + vi(t+Aty 



Xi(t + At) = Xi(t) + 



At, (57) 



where a,i(t) denotes the deterministic part of the right- 
hand side of Eq. (|3J), and z t ~ N(0,1) are independent 
realisations of a Gaussian distributed quantity with zero 
mean and unit variance. 

The velocity update ((56)) corresponds to decomposing 
the deterministic and stochastic parts of the accelerations. 
While the deterministic part corresponds to an Euler up- 
date, the stochastic part is a result of explicitely solving 
the stochastic differential equation 



dv 



(s) 



df 



6(0 



for the initial conditions v\ s \t) — at time t' = t. The so- 
lutions v^it + At) = VDAt zt are realisations of random- 
walk trajectories, which are Gaussian distributed with ex- 
pectation value zero, and variance DAt. 

Notice that, for sufficiently small update times, this 
update scheme should converge (in the statistical sense) 
to the true solutions of |2} and ([3]). In the simulations, 
we have set At = 0.04 s. To verify the convergence, we 
have also run some simulations with lower values of the 
time step (down to At = 0.001 s) and found less than 1% 
deviation. 



4.3 Gap distribution 

Figure [T] shows the simulated gap distributions for the 
stochastic OVM for two densities and two values of the 
symmetry parameter 7. For comparatively low densities 
(Fig. Q] (a)), both the predicted and the observed distri- 
butions are markedly asymmetric. Furthermore, the direc- 
tion of interaction plays a role as well. For the limiting case 
of a car-following model (7 = 0), the gap distribution is 
wider than in the symmetric (momentum-conserving) case 




60 80 100 120 140 160 180 200 

Gap (m) 




32 34 
Gap (m) 

Fig. 1. Stationary gap distributions for the stochastic optimal 
velocity model (sOVM) with the parameters specified in Ta- 
ble [l] on a ring road and densities of (a) p = 12veh.km, and 
(b) p — 30veh.km. Plotted are the simulated data (symbols), 
the theoretical distribution (|34|l (thick solid line), and, for the 
higher density, the Gaussian approximations (|42p (thin lines). 



7 = 1. Generally, there is a good agreement between the 
theoretical expressions and the data. The only exception is 
the large-gap tail for symmetric forces at the lower density. 
In contrast, for the car-following case 7 = 0, even the tails 
are reproduced correctly (within statistical fluctuations). 
The same aggreement has been found for the higher den- 
sity, irrespective of the value of 7. This is remarkable since 
additional assumptions have been necessary in Section [2] 
to derive the theoretical distributions for the car-following 
case. Consequently, one would expect larger errors com- 
pared to the isotropic case. 

Now we investigate the influence of the densities on 
the form of the gap distribution. Comparing Fig. [TJa) with 
Fig.rjjb), one may conclude that, when increasing the den- 
sity, the distributions become more and more symmetric. 
Further simulations showed that the distribution becomes 
significantly asymmetric if the single-particle kinetic en- 
ergy T = Dt/4 exceeds the effective potential energy 
[/(1/p) by at last one order of magnitude. Specifically, for 
the situation of Fig.[TJ we have T = 1 m 2 /s 2 for all values 
of p and 7 while the effective potential energy correspond- 
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Fig. 2. Stationary gap distributions for the stochastic power- 
law model (sPLM) with the parameter and simulation settings 
specified in Table [5] 



ing to plot (a) is given by U s ovm — (1 + 7) 0.127 m 2 /s 2 
and that of plot (b) by U s0 vm = (1 + 7) 95.0 m 2 /s 2 . 

For sufficiently high densities when the standard de- 
viation of the gap distribution is much smaller than the 
average gap 1/p, the Gaussian assumption (T4"2"l) should be- 
come valid. To determine the range of validity, we plotted 
the Gaussian approximation in the relevant Figure [ljb), 
in addition to the general distribution For the case 
7 = 1 corresponding to er 2 = 1.21m 2 , the Gaussian ap- 
proximation agrees nearly perfectly with the full theoret- 
ical curve (the curves overlap with no visible difference). 
For 7 = (<7 2 = 2.42 m 2 ), however, a significant difference 
is found, but the distribution already displays a sig- 
nificant skewness for this case. Further simulations showed 
that the Gaussian approximation is applicable whenever 
the full distribution (|34[) is sufficiently symmetric. 

In order to evaluate the robustness of the predictions 
with respect to different model types, we have simulated 
the gap distributions for the stochastic power-law model 
as well. The results are shown in Fig. [2] Apart from minor 
deviations at the large-gap tails, we found a remarkable 
agreement between theory and simulation. Moreover, the 
predicted distributions for the car-following case 7 = 
and the symmetric case 7 = 1 are significantly different 
both with respect to variance and shape, which justifies 
the particular specification ([9|) of the effective potential. 

Finally, we looked closer at the relaxation dynamics of 
the initially r)-correlated distributions, Eq. (|55p. towards 
the stationary distributions. A very slow relaxation could 
be a possible reason for the deviations found sometimes 
at the large-distance tails of the gap distributions. In Fig. 
we display snapshots of the evolution of the distribu- 
tion for different simulation times. The results show that 
the relaxation time is considerable, particularly for low 
densities. Moreover, the relaxation process is particularly 
slow at the tails, so it may be a plausible reason for the 
remaining deviations. 




60 80 100 120 140 160 180 
Scaled velocity deviation 



Fig. 3. Relaxation of the (initially 5-shaped) gap distribu- 
tion to its stationary distribution for the stochastic power-law 
model with the parameters specified in Tabled The figure com- 
pares simulation results at different points in time (symbols) 
with the theoretical stationary distribution (solid line). 
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Fig. 4. Stationary velocity distribution for the stochastic OVM 
with the parameters specified in Table [1] 

4.4 Velocity distribution 

In contrast to the gap distributions, the predicted veloc- 
ity distributions are always Gaussian. Moreover, the ve- 
locity variance should satisfy the fluctuation-dissipation 
theorem (|19[) . As a consequence, the variance may neither 
depend on the density nor on the direction of the interact- 
ing forces. Figure [4] shows simulated velocity distributions 
for the stochastic OVM at the higher density correspond- 
ing to Fig.QJb). One observes that, with the exception of 
small but systematic deviations from the Gaussian shape 
for the car- following case 7 = 1, all theoretical predic- 
tions are fulfilled. For the lower density p = 12 /km (not 
shown), the agreement was nearly perfect for all values of 
7- 

In contrast to the gap distributions, the agreement of 
the velocity distributions improves when going from the 
car-following to the conservative case and when decreas- 
ing the density. This can possibly be explained by the 
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Fig. 5. Stationary velocity distributions for the stochastic 
power- law model with the parameters specified in Table [2] The 
velocity distributions have been normalized with respect to the 
theoretical expectation value V , see Eq. (f54)l , and the variance 
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Fig. 6. Relaxation of the velocity distribution to the Gaussian 
stationary distribution for the stochastic power-law model with 
the parameters specified in Table [2] Initially, all particles had 
the same gaps s; = 1/p, and all velocities were equal to the 
expectation value. The figure compares simulation results at 
different points in time (symbols) with the theoretical station- 
ary distribution (solid line). 



distance from the instability point, see Eq. (fT3|) . For the 
densities p = 12 /km and 30 /km, the dimensionless dis- 
tances r — t c /t from the instability point are given by 
r = 0.001 and 0.132, respectively, while we have r = for 
the conservative case. Obviously, the agreement increases 
with the degree to which the requirement r <C 1 is satis- 
fied. 

In Figure we have plotted the simulated distribu- 
tions for the stochastic power-law model for different val- 
ues of the density and the directional parameter 7. In 
each simulation, we have normalized the distribution to 
the theoretical expectation value (fS"4")) and variance (fTT)|) , 
so we expect that all curves collapse onto each other in 



the ideal case. This collapse is, in fact, observed, thanks 
to values of r below 0.12 for all cases. 

Finally, we investigate the relaxation process from the 
(^-correlated intial velocity distribution to the stationary 
distribution. Figure[6]shows that there is a significant scale 
separation for the relaxation times: While the typical ve- 
locity relaxation time scale is of the order of seconds, it is 
of the order of hundred seconds for the gap distribution. 
After 1 000 s, even the tails of the velocity distribution are 
perfectly equilibrated, while for the gaps, this takes longer 
by a factor of more than one hundred. 

We conclude that, in contrast to the case of gap dis- 
tributions, long relaxation times cannot explain possible 
differences between the theoretical and simulated veloc- 
ity distributions, as noticeable in Fig. 2J These will be 
explained in the following. 



4.5 Kinetic energy and correlations 

One of the crucial assumption in the derivation of the 
gap and velocity distributions of Sec. [2] is the assump- 
tion of zero correlations, which requires that the system 
is far from any instability, i.e. r < 1. In classical ther- 
modynamic systems, it is well known and theoretically 
understood [16] that the energy contained in the fluctua- 
tions increases near a phase transition resulting in "critical 
opalescence" and other observable phenomena. The same 
has been found in driven thermodynamic systems such as 
Rayleigh-Benard convection or electrohydrodynamic con- 
vection [17j below the deterministic threshold, which will 
be further discussed in Sec. [3 

It is therfore very interesting to investigate the stochas- 
tic properties of our driven particle-systems as a function 
of the distance from threshold, i.e., varying the relaxation 
time from r = to r = r c or, cquivalently, the control 
parameter from r — to r = 1. 

Figure [7] shows the single-particle kinetic energy of the 
fluctuations as a function of the relaxation time for sev- 
eral values of the directional parameter 7. The kinetic 
energy has been normalized to the value resulting 
from the fluctuation-dissipation theorem ()19j) . As in the 
physical systems mentioned above, we found significantly 
increased, so-called "critical fluctuations" near the linear 
threshold, which is located at r c = 1.51 s for 7 = and at 
t c = 2.83 s for 7 = 0.2 while no such threshold exists for 
7=1. 

Finally, we compared the observed increase factor of 
the fluctation energy with the function l/yl — r of the 
scaled distance to the threshold (scaled control parame- 
ter). The agreement was astonishing for all investigated 
values of r and 7. 



5 Conclusion 

In this contribution, we have investigated the statistical 
properties of onedimensional dissipative driven many-par- 
ticle systems violating the law "actio = reactio". Such 



10 



Helbing and Treiber: Hamiltonian Statistics in Onedimensional Driven Dissipative Systems 



■B 2 



sOVM 






p = 30 /km 






y=0.0 • 1 






- approximation I 






y=0.2 a r 






approximation 1 






7=1.0 4 






equilibrium / 



















0.5 1 1.5 2 2.5 3 

Relaxation time x (s) 



Fig. 7. Average kinetic energy per particle for different val- 
ues of the relaxation time and three values of the anisotropy 
parameter 7. The amplification of the velocity variance with 
growing values of r reflects critical fluctuations close to the in- 
stability point, see Eq. (|13p . The energy has been plotted rela- 
tive to the theoretical value corresponding to the fluctuation- 
dissipation theorem (symbols). The solid curves display the 
function l/y/l — r of the reduced control parameter (|45[) . and 
the thin vertical lines give the asympotics for t-»t c (r-tl). 



systems can represent, for example, vehicular traffic or 
queuing systems with interactions. 

In the theoretical derivation, we have shown that such 
systems approximately show a Hamiltonian statistic when 
inter-particle correlations play no significant role. The the- 
oretical predictions were confirmed by simulations: With- 
out a single free parameter to fit, we quantitatively ob- 
tained the typical characteristic properties of Hamiltonian 
systems such as velocity and gap statistics corresponding 
to a canonical ensemble cx e~ n l e when the Hamiltonian Ti. 
contains the usual contributions of kinetic and potential 
energy as in physical Hamiltonian systems. Furthermore, 
the velocity variance satisfied the fluctuation-dissipation 
theorem. As only prerequisite, we have found that the sys- 
tem must be far away from any instability point. This is 
consistent with the theoretical requirement of vanishing 
correlations, since collective instabilities such as stop-and- 
go traffic correspond to highly correlated particles. 

In the first moment, these results appear to be quite 
surprising. For example, in traffic systems neither energy 
nor momentum are conserved during vehicle interactions, 
and the driving force keeps the system permanently far 
from equilibrium. While conservative systems conserve mo- 
mentum and energy in each single interaction, in the driven 
dissipative systems studied by us, the additional relax- 
ation term (v — Vi)/r causes the average velocity V to 
relax to the "free" or "desired" speed vq. 

However, even systems violating the law "actio=reactio' 
may be mapped to effectively conservative systems by a 
Galilei transformation. Defining motions relative to the 
stationary velocity, 

Ui=Vi-V (58) 



Eq. ^ becomes 

fill ' It ■ 

= ~ + /(*) - m/p) + 7 If(si-i) - f{l/p)] ■ (59) 

One sees that the constant terms —f(l/p) and 7/(l/p) re- 
sulting from the driving force and the relaxation dynamics 
supplement the interaction forces by counteracting forces 
- irrespective of the value of 7 - such as in momentum- 
conserving systems. 

One big difference of system not conserving momen- 
tum, however, remains when compared to conservative 
system: The conservative many-particle system always be- 
haves dynamically stable, while the dissipative system po- 
tentially produces stop-and-go waves, when the linear sta- 
bility condition [13] is not fulfilled. According to computer 
simulations, close to the instability point, the driven dissi- 
pative many-particle system tends to produce correlations 
between distances and velocities and between successive 
particles |18j . 

This corresponds to pattern formation phenomena that 
would not occur in conservative systems. Such kinds of 
pattern formation phenomena have, for example, been in- 
vestigated in fluid systems driven by thermal gradients 
(Rayleigh-Benard convection), coriolis forces (Taylor-Cou- 
ette flow), electrical fields (electroconvection), or concen- 
tration gradients (binary-mixture convection), see Ref. [17] 
for a review. Moreover, the increase of thermal fluctations 
when approaching a linear stability threshold from below 
has been investigated theoretically and experimentally for 
the above systems [19|20|21|22] . Near the threshold but 
in the regime of linear response, the fluctuations should 
increase according to a power law, where the scaling expo- 
nents depend on the dimensionality and symmetry classes 
of the systems [23]. Specifically, if the fluid systems are 
quasi-onedimensional, their fluctuations typically increase 
proportional to (1 — r) -1 / 2 where the reduced control pa- 
rameter r is defined in analogy to Eq. (|45| . with r replaced 
by a suitable driving force such as voltage or tempera- 
ture gradient. Very near the threshold, however, devia- 
tions have been observed empirically [24] , 

It appears that in our case the fluctuations scale pro- 
portionally to 1/^/1 — r as well. A theoretical founda- 
tion and a closer investigation of the above Fokker-Planck 
equation near the instability point and beyond will be 
subject of our future studies. 
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